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We compute f-mode travel-time sensitivity kernels for flows. Using a two-dimensional model, the scattered 
wavefield is calculated in the first Born approximation. We test the correctness of the kernels by comparing an 
exact solution (constant flow), a solution linearized in the flow, and the total integral of the kernel. In practice, 
the linear approximation is acceptable for flows as large as about 400 m/s. 
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I. INTRODUCTION 

The solar f modes are useful tools for studying perturba- 
tions that reside within the top 2 Mm below the surface of the 
Sun, and have been used as diagnostics as such before JU. 
The linear forward problem of time-distance helioseismology 
is to calculate functions, or kernels, which give the sensitivity 
of travel-time measurements to small-amplitude perturbations 
in the Sun. For a discussion of this general procedure, which 
is the one that we will follow here, see Gizon & Birch H. 

In this study we are interested in flow kernels, i.e. func- 
tions which give the sensitivity of travel-time measurements to 
steady, small amplitude, two-dimensional, spatially-varying 
horizontal flows. By flow kernel, we mean the function K 
defined by 



Tdiff(a;i|a;2) 



// 



u{x) ■ K(A,x)d 2 x, 



(1) 



where ram is the difference in travel time for waves propagat- 
ing in opposite directions between two fixed locations X\ and 
X2 on the surface, u — (u x , u y ) is the two-dimensional flow 
vector, and the integral is over the position vector x = (x, y) 
that spans the whole solar surface. We adopt a cartesian geom- 
etry for the sake of simplicity. The kernel K = (K x , K y ) is a 
two-dimensional vector which gives the sensitivity to flows in 
the x and y directions. Furthermore, we define A = X2 — x\ 
and the distance A = || A|| between the two points X\ and 
x 2 . 

In Section [IT] we describe our simplified model of f-mode 
propagation through an inhomogeneous moving medium. In 
Section [ill] we consider a uniform flow and compare the ex- 
act solution, its linearization in the flow, and the third-order 
approximation. In Section HVl we show example kernels. We 
discuss the results in the last section. 



II. ZERO-ORDER WAVEFIELD AND FIRST BORN 
APPROXIMATION 

In general, observations of solar oscillations are described 
by the filtered line-of-sight projection of the observed Doppler 
velocity v given by 



where the operator T describes the convolution with the point- 
spread function of the telescope and any additional filtering 
used in the data analysis, £ denotes the line-of-sight unit vec- 
tor, and cf> is the line-of-sight Dopplergram. 

As was done by Gizon & Birch |2|], we consider a constant 
density half-space with a free surface at z = 0. Gravitational 
acceleration, — gz, is constant. The medium is permeated by a 
steady, horizontal, divergenceless, and irrotational flow, u(x). 
The condition V m = is imposed by the continuity equation 
in the steady state. In addition, we require that V x u = to 
keep the problem two-dimensional. Under these conditions of 
smoothness, the fluid velocity remains irrotational at all times. 

As is appropriate for solar waves, we consider small- 
amplitude (linear) waves. At the free surface z = 0, the lin- 
earized dynamical and kinematic boundary conditions are 

p/p - 9V = n /p z = Q, (3) 
-iuri + u-Vr) = d z Q z = 0, (4) 

where is the wave velocity potential v = V0, f) is the ele- 
vation of the fluid's surface, p is a pressure perturbation, p is 
the background density, and IT is a stochastic pressure source 
that generates the waves. The above equations and the follow- 
ing ones apply to quantities that were Fourier-transformed in 
time (co is the angular frequency). In the bulk, the linearized 
momentum and continuity equations are 



iw6 



• ve 
v 2 e 



-p/p-T[Q], 



0, 



(5) 
(6) 



where T is a phenomenological frequency-dependent damp- 
ing operator Q2LI • 

In order to obtain the scattered wavefield computed in the 
first Born approximation, we expand any wave quantity q into 
a zero-order term q° (in the absence of the flow) and a first- 
order term 5q caused by the flow perturbation. In the bulk, 
the velocity potentials 0° and (50 both satisfy the three di- 
mensional Laplace's equation. Eliminating pressure and sur- 
face elevation, the zero-order and first-order surface boundary 
conditions at z = reduce to 



(d z - K )©° = = s°, 

99 



(d z - k)SQ = iu ■ V [(d u K)Q° + d^S ] = SS, 



ip(x,t) =T[t- 



(2) 



(7) 
(8) 

where k(u>) = to 2 / g + vyu)/g is the complex wavenumber 
at resonance and 7(0;) is a realistic damping rate [2]. In do- 
ing so, we have used the approximations <5IT = iu, • V<9 W IT 



2 



and (5r8° = i(d u ^y)u ■ V0°, which describe the perturba- 
tions to the source of excitation and the wave damping respec- 
tively (this can be understood by transforming to a co-moving 
frame). 

The solution of the zero- and first-order problems can be 
solved by introducing a Green's function. The final result for 
the observable ip (on the surface) is 

i/>°(x,u)=F[(£ b >V h + t t K)e ], (9) 
Sif>(x, uj)=T [(4 • V h + £ z k)5Q + i£ z (d w K)u • V6°] , 

(10) 

where i = (£h,l z ), Vh is the horizontal gradient, and the 
zero-order and first-order surface velocity potentials are given 
by 

e°(aj,cj) = 2vr J J G{x - x', u)S°(x', w)dV, (11) 
6Q(x, uj) = 2tt // G(x - x', uj)SS(x', uj)d 2 x', (12) 



with 



G(x,u) 



1 



(2tt) s JJ k-n{u) 



d 2 k, (13) 



where k = \\k\\ is the horizontal wavenumber. With the above 
expressions for the observed wavefield, one can compute the 
temporal cross-covariance function. 



III. UNIFORM FLOW: EXACT, FIRST-, AND 
THIRD-ORDER SOLUTIONS 

Before using the first-order Born approximation from the 
previous section to derive the sensitivity kernels (see Sec.lIVIi. 
we wish to consider here the case of constant flows, 



u = const. 



(14) 



In this case, <f>(x, t) can be calculated simply by considering 
the Galilean transformation x — > x — ut: 



(x,t) = (f>°{x - ut,*). 



(15) 



This transformation applies to all other quantitites in the prob- 
lem except for the filter T. So the filtered wavefield in Fourier 
space is then given by 



ip(k, uj) = F{k, w)</>°(fc, u-k-u), 



(16) 



where multiplication by F(k, to) accounts for the action of the 
filter operator T in Fourier space. The filter contains a mode- 
mass correction which is the ratio of the mode inertia in our 
model to the mode inertia in a standard stratified model. In 
general it also contains a phase-speed filter, but in the partic- 
ular case presented here where only f modes are excited, a 
phase-speed filter is unnecessary and therefore F is indepen- 
dent of u>. 



The expectation value of the cross-covariance function is 
given by |0] 

C{ Xi \x 2 ,uj) = J J P(fc,w)e ifeA d 2 fc, (17) 



where 



P(k,u) = E[\ip(k,cv)\ 



(18) 



denotes the expectation value of the power spectrum. In 
particular, the zero-order, or unperturbed (u = 0), cross- 
covariance is given by 



C ( Xi \x 2 ,cj) 



P v (k,u;)e ikA d 2 k. 



(19) 



Following Gizon & Birch [2], the travel-time difference can 
be measured from the cross-covariance function according to 



T^X^X-z) 



-2Rc J °° kjC°(A, to) AC(A, w) dw 
J™L0' 2 \C°{A,u;>)\ 2 du;' 



(20) 

using AC = C — C°, where Re takes the real part of the 
expression. The travel-time shift Tdiff (xi\x2), i.e., the differ- 
ence in time it takes for waves going from X\ to x 2 and for 
those going from x 2 to X\. In our problem, r^is is due to the 
flow. 

On one hand, the travel-time difference may be computed 
exactly for any value of u. This is simply done numerically 
by using eqs. (fT5T> - (f2Qb . 

On the other hand, an approximation to the travel-time dif- 
ference can be obtained by expanding the power spectrum as 
a Taylor series in u. This will provide us with a means to test 
the first-order Born approximation developed in Section[lll as 
well as to quantify higher-order terms. For a constant flow, 
the Taylor expansion of the power spectrum (eq. fT8l ) is 

P(k,co) =P°(k,uj)-k-u\F\ 2 d UJ \cl ) \ 2 
+ i(fc.«) 2 |F|Wl 2 
-i(fc-u) 3 |F|Wl 2 + --- 



(21) 



We may truncate this expansion after each successive term 
to study the first-, second-, and third-order travel-time differ- 
ences (using eq. l20l ). The second-order term does not af- 
fect travel-time differences because it only introduces pertur- 
bations to the cross-covariance that are symmetric in the time 
lag. 

Figure Q] shows the exact and approximate travel-time dif- 
ferences Tdiff /u versus A computed at fixed u = uA for 
u = 200 m/s and u = 400 m/s. We find that the first-order ap- 
proximation is within 10% of the exact value for A < 20 Mm, 
while the error is only about 1 % when terms up to u 3 are kept. 
For a flow that is twice as large, u = 400 m/s, non-linear ef- 
fects are more important. The non-linearity of Tdiff with u also 
increases as the separation distance A increases. 

Figure [2] shows the variations of Tdiff as a function of u 
at fixed distances A = 5 Mm and A = 10 Mm. The first- 
order approximation is reasonable for flows with amplitudes 
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FIG. 1 : F-mode travel-time differences versus distance for uniform 
horizontal flows. Two values are chosen for the flow amplitude, u = 
200 m/s and u = 400 m/s. The thick black solid line is the solution 
to the first-order approximation, the thin solid lines are the third- 
order solutions, and the dashed lines are the exact solutions. The 
green dots denote the values of the spatially-integrated kernels from 
Sec.HVl 



FIG. 2: F-mode travel-time differences versus uniform horizontal 
flow for different distances. Two distances are studied, A = 5 Mm 
and A = 10 Mm, given by the red and blue sets of curves respec- 
tively. For each set of curves, the thick solid line denotes the first- 
order solution, the thin solid line is the third-order solution, and the 
dashed line is the exact solution. 



less than about 400 m/s, although it gets worse for larger dis- 
tances. In particular this means that the first-order approxima- 
tion (and thus the kernels given in Sec. lIVI l is appropriate to 
study supergranular flows. 

IV. SENSITIVITY KERNELS 

We now return to the computation of sensitivity kernels in 
the first Born approximation in order to connect T^m with a 
spatially varying flow u(x). 

The first-order perturbation to the travel time is obtained by 
approximating AC in eq. d20j l by 

AC(ei,e 2 ,w) ~E [6ip*(x 1 ,u})ip°{x 2 ,Lu) 

+ r/>°*{xi,u)5^(x 3 ,u)]. (22) 

The linear dependence of AC and rais on u then follows from 
the expression of Sip (eq. |[T0l ). 

The kernel functions are identified according to their defi- 
nition, eq. ((T). The explicit expression of the flow kernels will 
be given in an upcoming publication. 

A pair of sensitivity kernels is given in Fig. [3] for a dis- 
tance A = 10 Mm (the observation points have coordinates 
xi = (5,0) Mm and x 2 = (-5,0) Mm). The kernel on 
the left, K x , gives the sensitivity to u x and the kernel on the 
right, K y , gives the sensitivity to u y . The flow kernel K x dis- 
plays elliptical and hyperbolic features, just as the kernels for 
the damping rate and the source strength derived earlier by 



Gizon & Birch [2]. The elliptical features have been called 
Fresnel zones in geophysics (e.g., JH]). The first Fresnel zone 
has a size (measured along x = 0) which is roughly given 
by (AA) 1 / 2 , where A ~ 5 Mm is the dominant f-mode wave- 
length. The hyperbolae are due to scattering of waves gen- 
erated by distant excitation events. For reasons of symmetry 
(the K y kernel is antisymmetric with respect to the lines x = 
and y = 0) the total integral of K y is zero. 

For consistency, we have checked that the total integral of 
K x is nearly exactly the same (within 0. 1 %) as the first-order 
approximation of T^g ju calculated under the assumption of 
a constant flow (Sec.lIiTland Fig . [TJ . This gives us confidence 
in the numerics. 

The effect of the line-of-sight projection of the Doppler ve- 
locity signal on the flow kernels are included in the model of 
Sec. [II] This is not presented here, although it is briefly dis- 
cussed in Jackiewicz et al. ||4J]. 



V. CONCLUSIONS 

Sensitivity kernels for flows have been calculated in the first 
Born approximation following the general recipe given by Gi- 
zon & Birch yfl. 

It has proven very useful to consider a uniform flow model, 
which allows one to find an exact solution for the travel times, 
as well as approximations to the exact solution in leading pow- 
ers of the flow u. 

It was shown that the first-order solution is equivalent to 
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FIG. 3: F-mode travel-time sensitivity kernels for flows. The separation distance of the observation points, denoted by the crosses, is A = 
10 Mm. The coordinates are x± — ( — 5, 0) Mm and X2 = (5, 0) Mm. The units of the colorscale are sMm~ 2 /(km/s). The K x and K y 
components of the kernel are shown in panel a) and panel b) respectively. 



the first Born approximation in the case of a uniform flow, 
and thus provides a test of the consistency of the kernels. This 
test revealed an accuracy within about 0. 1 %, giving us con- 
fidence in the numerical calculations. Furthermore, compar- 
isons of the first-order solution with the exact solution show 



that it is reasonable to study flows up to 400 m/s, with the 
caveat that this value depends on the distance A, with larger 
distances giving more error. Nonetheless, the linear model 
presented here is likely quite capable for the study of super- 
granular flows. 
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